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ABSTRACT 

We investigate the effects of the efficient production of cosmic rays on the evolution of supernova remnants 
(SNRs) in the adiabatic Sedov-Taylor phase. We model the SNR by coupling the hydrodynamic evolution 
with nonlinear diffusive shock acceleration (DSA), and track self-consistently the ionization state of the shock- 
heated plasma. Using a plasma emissivity code and the results of the model, we predict the thermal X-ray 
emission and combine it with the non-thermal component in order to obtain the complete spectrum in this 
energy range. Hence, we study how the interpretation of thermal X-ray observations is affected by the efficiency 
of the DSA process, and find that, compared to test particle cases, the efficient DSA example yields a smaller 
shock radius and speed, a larger compression ratio, and lower intensity X-ray thermal emission. We also find 
that a model where the shock is not assumed to produce CRs can fit the X-ray observational properties of an 
example with efficient particle acceleration, with a different set of input parameters, and in particular a much 
lower explosion energy. Additionally, we model the broadband non-thermal emission, and investigate what 
signatures result from the acceleration of particles. 

Subject headings: acceleration of particles - cosmic rays - shock waves - ISM:supernova remnants 



1. INTRODUCTION 

A supernova remnant (SNR) is the result of the interaction 
between supernova (SN) ejecta and the surrounding interstel- 
lar medium (ISM). As the SN ejecta expands it sweeps up the 
surrounding material, and forms a forward shock wave (FS) 
that compresses the ISM heating it to X-ray emitting temper- 
atures. The deceleration of the SN ejecta also produces a re- 
verse shock (RS) that travels inward from the contact inter- 
face, or contact discontinuity (CD), heating the ejected ma- 
terial. The evolution of SNRs and their emission character- 
istics can be divided into several different phases ( Chevalier 
|1977| l. Initially the expansion and thermal X-ray emission 
properties of an SNR are dominated by the ejecta. Once the 
majority of the ejecta energy has been transferred to the sur- 
rounding medium, the SNR enters the Sedov-Taylor phase. 
During this stage the hydrodynamic evolution can be approxi- 
mated by the self-similar adiabatic expansion of a shock wave 
originating at a poi nt explosion and pr opagating through an 
uniform medium ( |Sedov|[T959| |Taylor|[T950| l. For SNRs in 
this stage, the Sedov-Taylor similarity solution for the expan- 
sion of the blast wave, together with th e Rankine-Hugoniot 
relations (Rankine 1870 Hugoniot 1889), allow parameters 
of the SN and ambient medium to be estimated from thermal 
X-ray emission observations (Hamilto n et al.|1983| l. 

In addition to heating the surrounding medium, SNR shocks 
are thought to accelerate ambient particles to cosmic ray (CR) 
energies, through diffusive shock acceleration (DSA). Non- 
thermal X-ray emission has be en detected fr om young shell- 
type S NRs, including SN 1 006 (jKoyama et al.|1995)|Reynolds 

1998) , RX J1713.7-3 946 (jKoy ama et al.||1997[ |S lane et al. 

1999) , and Vela Jr. ( Aschenbach ll998| [siane et al.p OOl ). 
These X-rays are believed to be synchrotron radiation from 
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electrons accelerated to TeV energies at the SNR shock. Addi- 
tionally, if a significant fraction of the SN energy is placed in 
relativistic particles the hydrodynamic evolution of the SNR 
will be modified, and some evidence of this has been ob- 
served. X-ray observations of the relative positions of the FS 
and RS, and contact discontinuity (CD), in Tycho's SNR, and 
SN 1006 point to the shock compression ratios being modified 
by the acceleration of cosm ic-ray ions (Wa rren et al.||2003| 
Cassam-Chenai et al. 2008 ). Also, the postshock plasma tem- 
peratures observed in SNRs IE 0102.2-7219 and RCW 86, are 
lower than expected for their measured shock velocities, and 
hence als o indicate th at particl e acceleration a t their shocks is 
efficient ( |Hughes et al |2000| |HeIder et al.|2009| ). 

The evidence of efficient acceleration of cosmic rays at 
SNR shocks suggests that the DSA process must be incor- 
porated into the analysis of the evolution of remnants, and 
some early work has outlined the impact of particle ac celera- 
tion on observations of SN Rs in the adiabatic stage (Heavens 
1984||Dorfi & Bohringer|1993| >. Here, we investigate the ef- 
fects of efficient DSA on the evolution of SNRs in the Sedov- 
Taylor phase, and how it impacts the understanding of X-ray 
observations of remnants. To this end we use a hydrodynamic 
model coupled with nonlinear DSA, and with a self-consistent 
treatment of the non-equilibrium ionization (NEI) structure of 
the shock-heated plasma. 

In Section 3.2, we outline the Sedov-Taylor solution for 
the adiabatic expansion of SNRs, and the Rankine-Hugoniot 
equations for the propagation of shock waves. Additionally, 
we show how the relations between the physical parameters 
are modified when a fraction of the SN energy that goes into 
compressing the ambient medium changes, and the equation 
of state of the plasma is softened due to the acceleration of 
particles at the shock. We describe our model in more detail 
in Section 3.3, and present simulation examples of SNRs with 
both efficient and inefficient particle acceleration and discuss 
these results in Section 3.4. 

2. THE SEDOV-TAYLOR PHASE AND PARTICLE ACCELERATION 
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The Sedov-Taylor self-similar solution describes the expan- 
sion of a blast wave that originates in a point explosi on and 
propa gates in a homogeneous medium (Se dov||1959 Taylor 
|1950) . SNR shocks are believed to relax towards this descnp- 
tion once the mass of swept-up ISM is larger than that of the 
SN ejecta. In this section we present a simple analytical treat- 
ment that connects the SN explosion energy, the age, and the 
ISM density with parameters that can be obtained from ob- 
serving thermal X-ray emission from the shock-heated plasma 
in SNRs. 

For an SNR that transfers energy E into the heating of the 
surrounding plasma, with ISM pre-shock density po, the ex- 
pansion of the SNR at age t is described by the expression 



Rs 



a(j)Et 2 
Po 



1/5 



(1) 



where a is a function of the specific heat ratio of the gas 7, 
and is calculated numerically to solve the self-similar prob- 
lem ( |Sedov|1959| see Figure 75). However, as one expects a 
fraction 9 of the total SN explosion energy, Eg, to be placed 
in CRs, then we can define j3 = [a(l — 0)] 1//5 , and hence 
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From this result we can also derive the shock velocity to be 
v s = 2R s /(5t). 

The application of conservation laws across the shock front 
yield the Rankine-Hugoniot relations, which connect post- 
shock parameters, like temperature T and density p\, to the 
pre-shock conditions (McKee 1987, and references therein). 
For strong shocks, where the ratio of the shock velocity to the 
upstream signal velocity, or Mach number M, is much greater 
than 1, the compression ratio is given by 
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Furthermore, the pressure P$ of the post-shock material in this 
case is 
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where we can replace the pressure using the ideal gas law, 
P = pkT I (/xWh), where T is the temperature of the gas, and k 
is the Boltzmann constant, resulting in 
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In this expression p is the mean atomic weight, and me is the 
mass of the hydrogen atom. In the case of Sedov-Taylor SNRs 
we can express the shock velocity vs in terms of shock radius 
and age, and hence 
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In order to express the emission measure of the plasma 
£ = n e nnV, where n e and «h are the post-shock electron and 
hydrogen number densities respectively, in terms of the pre- 
shock conditions we need to consider the geometry of this 
simplified problem. The SNR shock in this basic picture is 



expanding spherically and compressing the ISM into a shell. 
From mass conservation considerations, the shell volume, V, 
must follow the relation p\V = po47r/?|/3, and using the com- 
pression ratio, as in Equation [3] 



V = 



(7-1)4 



nRi 



(7) 



(7+D 3 

Using Equations [3] and [7] we get the volume emission mea 
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where «o is the pre-shock proton number density. 

If the shock is not assumed to accelerate particles to CR 
energies then Equations [2] [6] and [HI can be used to estimate 
Eo, Po. and t, from the values of radius, plasma temperature, 
and emission measure, obtained from observations. In that 
case, when the fraction of the SN explosion energy that is 
placed in relativistic particles is 9 = 0, and the adiabatic index 
of the plasma is 7 = 5/3, and also assuming a mean atomic 
weight, p = 0.6, and a ratio of electron to hydrogen number 
density n c /«H = 1.23, we obtain the following expressions: 



no = 40.64 



■■ 0.075 



t = 423.9 



n_ 

10 60 cm" 



10 60 cm- 3 
kT s 



IkeV 



1/2 



1/2 



-1/2 



Rs_ 
lpc 

Rs_ 
lpc 



-3/2 



yr. 



kT s \ 
IkeV J 



Rs_\ 
lpc) 



3/2 



(9) 
(10) 

x 10 51 ergs . 



(ID 

Cosmic ray acceleration at the shock modifies this analy- 
sis, since some of the explosion energy goes into relativis- 
tic particles, and hence 6 > 0. Efficient particle acceleration 
also makes this a two-fluid problem, with a thermal plasma of 
adiabatic i ndex 7 = 5/3, an d a relativistic gas of CRs with 
7 = 4/3. I Chevalier] ( |1983| uses enthalpy arguments to ad- 
dress this problem analytically, and defines an effective spe- 
cific heat ratio, 

7*=^, (12) 
3(l + w) 

where w = p cr /p s is the fraction of the post-shock pressure in 
relativistic particles, which is dependent on the fraction of the 
total energy going into particle acceleration. Additionally, a 
fraction of the high energy particles accelerated at the shock 
is expected to escape from the shock, furt her softening the 
equat ion of state in the shocked plasma (Berezhko & Ellison 
1999). 

The softening of the equation of state of the plasma, and 
lower energy available for plasma heating, result in higher 
compression ratios and lower post-shock temperatures. Fig- 
ure 1 illustrates how the radius, the temperature and the emis- 
sion measure, calculated using Equations |2| |6| and [8] change 
as a function of the fraction of SN explosion energy placed 
into CRs. The variation of the effective adiabatic index of the 
post- shock plasma with 9 is determined using the results from 
|Chevalier| ( [T983) . 

In Figure |2J we show how the retrieved values of explo- 
sion energy, ambient density, and age of the SNR, obtained 
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9 = E cr IE 

Figure 1. Profile of the radius (solid), temperature (dotted), and emission 
measure (dashed) as functions of the energy placed in CRs, relative to the 
values of these parameters for the case where no efficient acceleration is tak- 
ing place (8 = and 7 = 5/3). 



using Equations |9|[TT| which assume that no CR acceleration 
has taken place, deviate from the input parameters used in 
Equations |2j [6] and [8] if the shock places some of its en- 
ergy into the acceleration of CRs, i.e., 6 > 0. This simple 
approach is a useful illustration of the approximate modified 
trends of the Sedov-Taylor analysis of the evolution of SNRs 
where efficient particle acceleration at the FS is taking place. 
However it does not include several important factors in the 
evolution of SNRs such as the free expansion phase, when the 
dynamics of the remnant are dominated by the expansion of 
the SN ejecta, and the temporal variation of the effective adia- 
batic index of the plasma. Nonetheless, the plot identifies the 
expected trends for parameters derived under the assumption 
that no CR acceleration has occurred. 
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Figure 2. Approximate trends of the inferred SN explosion energy Eq, age 
f, and ambient proton number density hq, as a function of the fraction of 
the SN explosion energy placed in relativistic particles, assuming no efficient 
acceleration is taking place (6 = and 7 = 5/3). 



3. THE CR-HYDRO MODEL WITH X-RAY EMISSION 

The CR-hydro model is a hydrodynamic simulation of an 
SNR that incorporates particle acceleration at the shock, and 
determines the ionization state of the shocked plasma at ev- 
ery time step. The X-ray emission profile is then calculated, 
using the NEI information for the thermal component, and 
the electron spectrum for the simulated s ynchrotron em ission. 
This model is described in more detail in Ellison et al. ( |2007| ), 
|Patnaude etaLl ( |2009] l; |Patnaude e"taLl ( |2010[ ), and references 
therein. 

CR-hydro calculates the hydrodynamic evolution of an 
SNR with a spherically symmetric model described in 
et al. (2007), and refe rences therein. A semi analy tica 
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of non linear DS A by Amato & Blasi ( 2005 ) and Blas i et al.| 
(2005), is coupled to the hydrodynamic calculation. Given 
an input DSA efficiency of the shock £dsa, the model deter- 
m ines the appropr iate injection parameter (£ in Equation [25] 
in |Blasi et al.|2005] l as the system evolves. The injection pa- 
rameter is then used to solve the DSA problem, together with 
the shock information and ambient density, temperature, and 
magnetic field. Hence, the energy loss from escaping accel- 
erated particles is removed from the shocked plasma, and the 
adiabatic index of the shocked gas, 7 s k, is determined from 
the particle distribution function, which includes the appro- 
priate ratio of relativistic and non-relativistic particles. This 
adiabatic index is then used in the hydrodynamic calculation, 
coupling DSA to the shock evolution. 

CR-hydro tracks the time-dependent hydrodynamics of the 
system in Lagrangian mass coordinates, subdividing it in a se- 
ries of concentric shells of shocked material. As the mass of 
each shell remains constant, the volume of the shell as the 
system evolves is adjusted to account for variations of the 
pressure. Both the RS and FS are included in the model, and 
shocked ejecta and shocked ISM are separated by a contact 
interface. As the shocks overtake new material, new shells 
are added to the system. A schematic view of the geometry of 
the model is shown in Figure [3] 

At each time step, the NET state of the shocked plasma, 
between the FS and CD, is calculated self-consistently, us- 
ing the ionization structure, free electron number dens ity, and 
electron temper ature. An updated version of the Raymond 
( |1977 1 plasma emissivity code is then used to ob 
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tain the thermal X-ray emission from the system (Brickhouse 
|et al.| [T995 ), where the resulting NEI fractions of heavy ele- 
ments, and plasma temperature an d de nsity, are used as inpu t 
parameters. Patnaude et aL|( |2009[ l, andlPatnaude et al. (2010), 
provide a detailed description of this NEI model. 

The resulting electron spectrum from the acceleration pro- 
cess is used to calculate the synchrotron radiation from the 
system, as well as the IC emission assuming standard seed 
photons from the cosmic microwave background. Addition- 
ally, the relativistic Bremsstrahlung photon spectrum is ob- 
tained using the spectral distribution of accelerated electrons. 
The model also calculates the pion decay emission spectrum, 
resulting from proton-proton collisions, based on the popula- 
tion of protons accelerated to relativisti c speeds at the FS, an d 
using the parameterization described in Ka mae et al.| f2006). 

An important factor determining the shape of the syn- 
chrotron spectrum is the nature of the magnetic field. Because 
the compression of the fluid is altered by DSA, the compres- 
sion of the magnetic field is also expected to be affected by 
the efficiency of this process. The ad hoc compression model 
of the magnetic field used in this simulation is described in 
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Figure 3. Diagram of the SNR model used in CR-hydro, which shows the 
positions of the FS, RS, and CD. The model is divided in concentric spherical 
shells of shocked material overtaken by the FS, which evolve through time. 



|Ellison & Cassam^Ch enai (2005). This treatment uses the far 
upstream magnetic field to calculate the Alfven heating rate 
in the shock precursor, which in turn results in changes in the 
subshock compression and the acceleration efficiency of the 
system. Magnetic field amplification (MFA) at the FS is now 
believed to be an importa nt effect in young S NR shocks (e.g., 
Vink & Laming 2003 Uchiyama et al.|2007), and is included 



in the model as an ad hoc magnification factor, set explicitly as 
an input parameter, B amp . The compressed field immediately 
behind the shock is increased by a factor B amp . The resulting 
compressed field is used to calculate the synchrotron emission 
from the accelerated electron population. 

We study a set of models where we explore different values 
of the ejecta kinetic energy from the SN explosion £sn, and 
the acceleration efficiency, £dsa- In all examples in this work 
the model is initialized using input parameters similar to those 
of Type la supernovae, i.e, mass of the e jecta M e \ = 1.4Mp , 
exponential density profile of the ejecta (Dwarkadas 2000), 
and kinetic energy in ejecta from the SN explosion Eq = 10 M 
erg. We assume a uniform ambient proton density n p , and 
magnetic field Bo = 5/iG. Additionally, we select explicitly 
the value of the relativistic electron to relativistic proton ratio 
K ep = 2%, in line with estimates inferred from direct mea- 
surements of cosmic ray spectra at Earth. This parameter de- 
termines how thermal and non-thermal X-ray emission scale 
relative to each other, as well as whether the 7-ray emission is 
dominated by leptonic or hadronic emission processes. Fur- 
ther discussion of the input p arameters require d fo r the CR- 
hydro mode l is provided in Ellison et al. (2007 ) and Pa tnaude| 
|et al.| ( [2009] ). In our next work, we intend to apply the model 
presented here to real cases of SNRs, for which it will be cru- 
cial to explore the parameter space extensively. 

The thermal X-ray emission model obtained from the 
plasma emissivity code is combined with the non-thermal 
emission in this range to obtain the complete X-ray emission 
profile. In order to simulate actual observations of SNRs in 
X-rays, the resulting spectrum model is then folded through 
the instrument response of the Advanced CCD Imaging Spec- 
trometer (ACIS) on board the Chandra X-ray Observatory . 
The distance to the SNR is fixed at 5 kpc, and the Galactic ab- 
sorbing column density is set to be A^h = 5 x 10 21 atoms cm 2 . 
No background contribution is included in the model or the 



subsequent analysis. 

4. RESULTS AND DISCUSSION 

In order to investigate how SNRs in the Sedov- Taylor phase 
are affected by efficient particle acceleration, we study a set of 
examples where the model is allowed to evolve to ages where 
the swept up ISM mass is much greater than the ejecta mass. 
First, we explore the test particle (TP) case, where the parti- 
cle acceleration process is not efficient; then, using the same 
set of input parameters, we study the case were DSA is effi- 
cient; and finally, we consider how the efficient case can be 
interpreted, and compare its results to those of test particles 
models. 

4.1. Test particle and efficient acceleration cases 

We consider the TP case by setting the DSA efficiency of 
the model to £dsa = 0.1%, and the ambient proton density to 
n p = 1 cm" 3 . The model is allowed to evolve to ?snr = 10,000 
yr, and at the end of the simulation we obtain the following 
parameters: FS radius /?fs = 12.6 pc; FS speed Vps = 510 km 
s , overall FS compression ratio R tot = 4.0; CD radius Rcd = 
8 .5 pc; and a swept up mass M sw = 280M Q . The fraction of the 
total explosion energy placed into cosmic rays in this example 
is 0=1.6%. 

The expected value of the radius of the FS for this set of 
parameters can be derived from Equation[T] which yields = 
12.5 pc. The small deviation from this value obtained from 
the CR-hydro model (>1%) can be attributed to the fact that 
the simple form of the Sedov-Taylor solution that results in 
Equation [T] does not take into consideration the initial phase 
of the SNR, when the SN ejecta freely expands into the ISM. 

For the efficient acceleration case we assume that 40% of 
the shock ram kinetic energy flux is placed in superthermal 
particles, i.e., £ D sa = 40%. At an age of 10,000 yr, and for 
« p = 1 cm -3 , the output parameters for this example are: FS 
radius /?ps = 12.1 pc; FS speed Vps = 480 km s" 1 ; overall FS 
compression ratio R tot = 4.9; CD radius /?fs = 8.6 pc; and a 
swept up mass M sw = 250M Q . In this case 28% of the 10 51 erg 
SN explosion energy goes into the acceleration of particles at 
the shock. 

The thermal X-ray emission is controlled in a complex way 
by the efficiency of the process of cosmic ray production. To 
understand this situation in more detail, the temporal evolu- 
tion of the radius, proton and electron temperatures, volume 
emission measure, and thermal X-ray flux (in the range 0.3- 
5.0 keV), for the TP and DSA efficient cases are presented in 
Figure |4] For the input parameters studied, the forward shock 
radius (top panel) and particle temperatures (second panel) of 
the TP case are consistently larger than those of the efficient 
acceleration case. While the emission measure is expected to 
broadly scale as the compression ratio times swept up mass, 
in a complex hydrodynamic model such as the one presented 
here, the value of the total emission measure will also de- 
pend on how the shocked gas evolves after the initial com- 
pression, and on how the compression ratio varie s with time. 
[ Ellison et aLl ( |2007| ), and |Patnaude et al.| ( |2009l i; |Patnaude et| 
|al.| ( |2010| i find that for ages t < 1Q 3 years, the volume emis- 
sion measure increases with DSA efficiency. The third panel 
in Figure [4] shows that, for the parameters used, the volume 
emission measure of the DSA efficient case is indeed larger 
than that of the TP case at 1 , 000 years, yet the ratio edsa/ctp 
decreases with age, and at approximately 4.000 years it be- 
comes smaller than 1. Figure [4] (bottom panel) shows that 
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the thermal X-ray flux from the TP case becomes increas- 
ingly more intense with age than that of the DSA efficient 
case, which is the combined result of the temporal evolution 
of the emission measures, electron temperatures, and ioniza- 
tion conditions. 
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Figure 4. The top panel shows the variation of the radius of the forward 
shock as a function of remnant age. As in all panels, the red solid line rep- 
resents the TP model, and the efficient acceleration case with £dsa = 40% is 
shown as a black dashed line. The emission measure weighted mean tem- 
peratures of protons and electrons are presented in the second panel. The 
third panel shows the evolution of the total volume emission measure of the 
shocked ISM, as well as the ratio of this parameter between the efficient and 
TP cases, presented as a blue dotted line (right axis). The bottom panel 
presents the thermal X-ray flux in the energy range 0.3-5.0 keV. The SNR 
models shown have input parameters Eq = 10 5 ' erg, and n v = 1 .0 cm -3 . 



Figure [5] shows the X-ray emission spectra generated from 
the TP example (top panel), and the efficient acceleration 
case (bottom panel). In both plots, the solid histogram rep- 
resents the thermal emission component, and the dotted line 
shows the contribution from non-thermal emission in this en- 
ergy range. The thermal emission spectra are obtained us- 
ing the plasma emissivity code with the NEI information ob- 
tained from CR-hydro for each example. Solar abundances 
are adopted for C, N, O, Ne, Mg, Si, S, Ar, Ca, and Fe, and 
the helium number density is fixed at 9.77% the proton num- 
ber density n p . 
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Figure 5. X-ray spectral energy distribution from the TP model (top panel), 
and efficient acceleration case with £qsa = 40% (bottom panel). The SNR 
models shown have input parameters t = 10 4 yr, Eq = 1.00 X 10 51 erg, and 
n p = 1.0 cm -3 . The adopted distance to the SNRs in both cases is 5 kpc. The 
solid line histogram represents the thermal component of the emission, and 
the dotted line histogram shows the predicted non-thermal emission in this 
energy band. 



The thermal component of the spectrum in the TP case 
clearly dominates over the non- thermal emission in the X-ray 
band. The total X-ray luminosity (in the 0.3-10 keV band) 
is Lx = 2.0 x 10 37 erg s , and the non-thermal contribution 
only represents ~ 10~ 3 % of this amount. Similarly, in the ef- 
ficient acceleration case thermal emission is also dominant, 
however the contribution from the non-thermal spectrum is 
much more significant that in the TP example, and at E > 5 
keV it overtakes the thermal component, apart from a few 
strong emission lines. The X-ray luminosity in this case is 
Lx = 1.2 x 10 37 erg s" 1 , and the non-thermal emission repre- 
sents a 0.9%. 

In Figure [6] we show the X-ray emission models folded 
through the Chandra ACIS instrument response. The sim- 
ulated X-ray spectrum, for a 5 ks observation of the entire 
SNR, from the TP case is shown as red crosses, and the ef- 
ficient acceleration case is presented as black crosses. Both 
spectra are clearly dominated by the thermal emission, and 
emission lines are very apparent. However, the emission from 
the TP case is obviously more intense than that from the case 
with efficient acceleration. 



4.2. Interpretation of the efficient acceleration case 

To investigate the scenario in which the emission from 
an SNR with efficient particle acceleration is treated with a 
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Figure 6. Simulated Chandra ACIS X-ray spectra of the model SNRs, with 
integrated 5 ks observation time. A Galactic column density A^h = 5 x 10 21 
cm~ 2 is used, and the adopted distance is 5 kpc. The SNR models shown 
have input parameters t = 10,000 yr, Eq = 10 5 ' erg, and n p = 1.0 cm -3 .The 
red crosses represent the TP case (edsa = 0. 1%) and the spectrum of the case 
with efficient DSA (edsa = 40%) is shown as black crosses. 



Sedov analysis that neglects the effects of acceleration, we 
searched for a TP case that fits the spectrum (as well as the 
radius) obtained in the efficient acceleration case by investi- 
gating variations in the input parameters such as age and SN 
explosion energy. The resulting X-ray models were then com- 
pared to the simulated Chandra spectrum resulting from the 
efficient acceleration case detailed in Section 4.1. 

An example of a TP case with results similar to those of the 
efficient acceleration case, is the model for a SNR with input 
parameters: s = 0.1%; age t = 12,600 yr; E Q = 0.56 x 10 51 
erg; and ambient proton number density n p = 1.1 cm" 3 . The 
resulting fraction of the explosion energy placed into CRs is 
9 = 1.8%, and the radius of the FS is R FS = 12.1 pc. 

In Figure U\ we show the model spectrum (red histogram) 
from this TP example, together with the simulated X-ray ob- 
servation of the efficient acceleration case, for a 5 ks exposure 
(black crosses). The best-fit histogram (shown in black) is ob- 
tained assuming the model to be an absorbed combination of 
the thermal emission model from the simulation and a small 
powerlaw component (green histogram). The resulting best- 
fit model is obtained with parameters Nu = 5.3 x 10 21 cm" 2 , 
powerlaw index T = 3.89, and unabsorbed powerlaw flux (in 
the 0.3-10 keV band) f x =4x 10" n erg cnrV 1 . The abun- 
dances of N (1.4 relative to solar), and S (1.2 relative to solar) 
were modified slightly to improve the fit. The reduced % 2 
statistic of the model is 1.27 for 147 degrees of freedom, and 
the statistical variations of the data from the model are shown 
in the lower panel in Figure [7] The TP fit appears quite rea- 
sonable, and this example illustrates that spectra from evolved 
SNRs where efficient CR acceleration took place can be well 
modeled by TP emission models, but the parameters derived 
from such fit would differ from the correct values. In this case 
the derived explosion energy is lower, and the derived density 
and age are higher than the input parameters, which is consis- 
tent with the trends illustrated in Figure [2] 

With the purpose of expanding the parameter space stud- 
ied, this analysis was repeated for a similar set of initial con- 
ditions, but considering expansion into an environment with 
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Figure 7. Simulated Chandra ACIS X-ray spectrum from the efficient accel- 
eration case (black crosses), with integrated 5 ks observation time, with input 
parameters as in Figures |jand|6] A Galactic column density Nu = 5 x 10 
cm~ 2 is used, and the adopted distance is 5 kpc. The best-fit model (black his- 
togram), is the absorbed (Nu = 5.3 X 10 21 cirT 2 ) combination of a powerlaw 
component (green histogram), with the emission model for an SNR model 
with input parameters t = 12, 600 yr, Eg = 0.56 X 10 51 erg, and rip = 1.1 cm" 3 
(red histogram). The bottom panel shows the statistical variations of the data 
in reference to the model. 



lower proton density, n p = 0.1 cm" . At 10,000 years, the for- 
ward shock of the DSA efficient case (e = 40%) expands to 
Rfs, = 19 pc, and the fraction of the total explosion energy de- 
posited into cosmic rays is 9 = 68%. We find that the X-ray 
emission from this model is satisfactorily fit by that of a test 
particle case, like in the models considered above. The TP 
model input parameters required to reproduce the DSA effi- 
cient case are e = 0.1%; age t = 12,000 yr; E = 0.54 x 10 51 
erg; and the same ambient proton number density n p = 0.1 
cm" 3 . This model also expands to Rps = 19 pc, and in Figure 
[8| we show the model spectrum (red histogram) from this TP 
example, together with the simulated X-ray observation of the 
efficient acceleration case, for a 20 ks exposure. The result- 
ing best-fit model is obtained with parameters Nu = 5.1 x 10 21 
cm" 2 , and no powerlaw component. The abundances of N and 
S were varied slightly relative to solar values, as in the higher 
density case. The reduced x 2 statistic of the model is 2.1 for 
148 degrees of freedom, and the statistical variations of the 
data from the model are shown in the lower panel in Figure [8] 

These results are reminiscent of the low SN explosion ener- 
gies (below than the canonical value of 10 51 erg) that have 
been estimated fro m X-ray observ ations for several SNRs, 
e.g., G272 , 2-3.2 ( |Harrus et al.pXJT) , and G299.2-2.9 ( |Park"| 
|et al.|2007j > . Our study suggests that these estimates, inferred 
from Sedov-Taylor analysis of remnants, could be the result 
of the modification of the SNR shock due to particle acceler- 
ation, and not underenergetic explosions. 

The broadband spectrum of the efficient acceleration exam- 
ple, with density n v = 1 cm" 3 , is shown in Figure m where 
we plot the thermal X-ray model, as well as each orthe non- 
thermal emission components (synchrotron, IC, non-thermal 
bremsstrahlung, and pion decay), and the summed non- 
thermal spectrum. Additionally, the summed non-thermal 
emission model from the test particle case modified to fit the 
X-ray spectrum of the efficient DSA case is also presented. 
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Figure 8. Simulated Chandra ACIS X-ray spectrum from an efficient ac- 
celeration case (black crosses), with integrated 20 ks observation time, with 
model input parameters / = 10,000 yr, Eo = 1.00 X 10 51 erg, and n p = 0.1 
cm~ 3 . A Galactic column density Nn = 5 X 10 2 ' cirT 2 is used, and the 
adopted distance is 5 kpc. The best-fit model (black histogram), is the ab- 
sorbed (Nn = 5.1 X 10 21 cm -2 ) emission model for an SNR model with in- 
put parameters t = 12,200 yr, Eo = 0.54 X 10 5 ' erg, and n p = 0.1 cm -3 (red 
histogram). The bottom panel shows the statistical variations of the data in 
reference to the model. 
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Figure 9. Broadband spectrum from the model with £dsa = 40%, no = 1 
cm -3 , £sn = 10 5 ' erg, and t =10,000 yr. The black histogram indicates the 
thermal X-ray emission, and the blue solid line represents the synchrotron 
component. Inverse Compton radiation is shown in green, pion decay in red, 
and the magenta curve represents the non-thermal bremsstrahlung emission. 
The dotted blue line is the sum of the non-thermal emission components for 
the DSA efficient case, and the cyan curve shows the sum of the non-thermal 
radiation from the TP case, where no = 1.1 cnr 3 , E$^ — 0.56 X 10 51 erg, 
and t =10,000 yr. The grey bowtie indicates the integrated sensitivity of the 
Fermr'-LAT after 1 year. 



While the thermal X-ray emission from the case with efficient 
cosmic ray acceleration, and that of the test particle case are 
similar (as shown in Figure [7), the non-thermal emission is 
clearly enhanced by the DSA efficiency. Figure |9] also shows, 
as a grey bowtie, the sensitivity after 1 year of observations of 
the Large Area Teles cope (LAT) on board the Fermi Gamma- 
ray Space Telescope ( [Atwood et al. [ 2009 ). Non-thermal emis- 
sion from the DSA efficient case would clearly be detectable 
in the MeV-GeV range with the LAT, while the TP case would 
not bright in that energy band. 

5. CONCLUSIONS 

We have presented a set of models of SNRs in the adiabatic 
phase, and considered the effects of efficient DSA on their 
hydrodynamics and emission characteristics. In this work we 
study only a limited range of parameters, yet it is clear that 
the production of CRs by the SNR shock significantly modi- 
fies the system. Compared to the TP case, the efficient DSA 
model yields smaller shock radius and speed, larger compres- 
sion ratio, and lower intensity X-ray thermal emission. We 
also found that a model where the shock is not assumed to 
produce CRs can fit the X-ray observational properties of an 
example with efficient particle acceleration, with a different 
set of input parameters, and in particular a much lower explo- 
sion energy. 

The analysis of observations of thermal X-ray emission 
from SNRs is often done using the Sedov-Taylor interpreta- 
tion, and assuming that shocks do not place a significant frac- 
tion of their energy into CRs. SN explosion energies lower 
than the canonical value of 10 51 erg have been estimated from 
X-ray observatio ns of several SNRs, including G272 .2-3.2 
( |Harrus et al.pjOl) , and G299.2-2.9 ( |Park et al.|2007) . Our 
study suggests that low SN explosion energies inferred from 
Sedov-Taylor analysis of remnants could be the result of the 
modification of the SNR shock due to particle acceleration. 

Compared to the TP case, non-thermal emission is much 
more intense with efficient acceleration (see Figure|9]l. There- 



fore, broadband modeling of SNR observations appears to be 
crucial for a appropriate understanding of the remnant's char- 
acteristics and those of its surroundings. 
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